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We demonstrate an experimental method for direct 2D and 3D imaging of magnetic radio-frequency (rf) field 
distribution in metal-dielectric structures based on traveling wave (TW) magnetic resonance imaging (MRI) at 
ultra-high field (>7T). The typical apparatus would include an ultra-high field whole body or small bore MRI 
scanner, waveguide elements filled with MRI active dielectrics with predefined electric and magnetic properties, 
and TW rf transmit-receive probes. We validated the technique by obtaining TW magnetic-resonance (MR) images 
of the magnetic field distribution of the rf modes of circular waveguide filled with deionized water in a 16.4 T 
small-bore MRI scanner and compared the MR images with numerical simulations. Our MRI technique opens up 
a practical non-perturbed way of imaging of previously inaccessible rf field distribution of modes inside of various 
shapes metal waveguides with inserted dielectric objects, including waveguide mode converters and transformers. 


I. INTRODUCTION 

Imaging of the rf magnetic field distribution is important 
for various applications that involve designing waveguide el¬ 
ements, mode splitters, and mode converters, where specific 
rf (microwave) field profile is desirable for engineering high- 
power rf propagation over the long distances, microwave 
heating in various research and commercial applications “ . 
Most of the measurements of the static magnetic field are 
done by utilizing various probes and consecutive mapping 
the field point by point. Electromagnetic field patterns cur¬ 
rently can be imaged from surface radiometry by using ther¬ 
mographic pattern analyzers ■ . Previously the dielectric 
waveguides were studied in optical band and corresponding 
electromagnetic mode patterns were observed in the optical 
fibers ■ . However modes imaging in rf band is still tech¬ 
nologically challenging. Recent advents of ultra-high field 
MRI make it possible to propagate resonant rf waves in over- 
moded waveguide regime inside MRI scanners . The prop¬ 
agating rf transverse magnetic field component becomes a 
spin excitation field for an MRI active medium thus imprint¬ 
ing a map of the transverse Hi field distribution (in MRI a 
flux density B\ is commonly used instead of H-field) into a 
sample. Such an image of the field map can be obtained 
in real time provided that the waveguide is physically com¬ 
patible with the high DC magnetic field inside the bore of 
MRI scanner. A complimentary way to obtain an rf field 
distribution is by computer simulation using various electro¬ 
magnetic field solvers that require a lot of computational 
resources and multiple solvers for solution validation. In our 
work we developed and demonstrated a TW MRI technique 
to directly image rf modes inside a circular waveguide filled 
with MR active high-permittivity dielectrics. 

In MRI, Hi field maps of near-field coils at the reso¬ 
nance frequency of the MR scanner are obtained using a 
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uniform dielectric sample in a process called rf shimming 
for optimization of custom-built coils ■ . Therefore, the 
modes of such volume near-field coils are well character¬ 
ized through both computer simulations and MRI. In recent 
work, MRI images of the modes of cylindrical dielectric res¬ 
onators placed in high-field MR scanner were reported ■ 
Although, conceptually similar to our method the latter re¬ 
lies on dielectric resonances of isolated cylinders with exact 
dimensions that correspond to specific dielectric resonances. 

At high field MRI (> 4 T) the propagating waves ef¬ 
fects inside dielectrics become important 14-16 and can be 
used for imaging provided that the cut-off requirements for 
the mode propagation is fulfilled. Traveling wave (TW) 
MRI is a far-field imaging technique that relies on rf mode 
propagation in a cylindrical magnet bore utilized as an rf 
waveguide. This approach was originally experimentally 
demonstrated at 7 T and at clinically relevant field of 
3 T 1 . However, the main challenge for the TW propaga¬ 
tion in most conventional MRI scanners and spectrometers 
is the relatively small bore radius-to-wavelength ratio (a/A) 
even at ultra-high field strengths. The cut-off condition in 
the cylind rical bore is giv en by the dispersion relationship: 
h mn = yW' - (x^ n /a) 2 , where h is the propagating vec- 
tor, m, n are integer numbers that are related to azimuthal 
and radial field variations, respectively, k is the k-vector of 
free propagating wave inside a dielectric (with s,p), a is 
the radius of the guide, x'mn ' s the n- th root of the deriva¬ 
tive of the Bessel function J m of the first order m (m=0, 
1, 2). So the cut-off frequency for the lowest TE mode is 
f c = x^ nn /(27ra- v //i5) and the critical wavelength of a mode 
is A c = c// c . For example, in a 7 T human-size MRI system 
with a metal guide fitted to the bore (maximum diameter 
of 60 cm) at a resonant frequency of 298 MHz (free space 
wavelength A = 100.6 cm), only a single mode ( TEn ) will 
propagate in an unloaded bore. Although most of TW MRI 
were done in a large bore MRI scanner at 7 T, rf imag¬ 
ing might be impractical there due to relatively small field 
of view (FOV ^ 50cm) to the rf wavelength ratio ( e.g ., 
A = 30 cm inside a dielectric with £ r =80). At high polariz- 
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FIG. 1 . Setup: a) an MRI scanner comprised by the metal 
bore with an image encoding gradient insert, a dielectric guide 
(DG): acrylic tube filled with Di-water, FOV=held-of-view, 
and b) a loop-coil (LC) oriented in axial plane. The rf mode 
propagates several wavelengths to the imaging FOV (dashed 
rectangle); its location with respect to the end of the tube 
can be varied. 


ing Bq fields (high resonance frequency) the rf wavelength 
inside a waveguide with a dielectric insert that has a high 
e r becomes small enough to facilitate propagating modes 
of such a hybrid waveguide. In our work we incorporated 
a dielectric waveguide that fits into FOV of a 16.4 T (698 
MHz) MRI scanner with a metal screen formed by the gradi¬ 
ent insert itself and devised the transverse rf magnetic field 
distribution of such waveguide from the intensity pattern of 
MR images. 


II. METHOD AND THEORY 


The complete solutions for the components of H field for 
TE (TM) waves inside a hollow circular waveguide are given 
by analytical expressions below 
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where c m is a constant. The field for TM modes can be de¬ 
rived from the Helmholtz equation by setting H z = 0 and 
the dispersion relationship is modified by substituting x'mn 
for Xmn ~ the n-th zero of the Bessel function. Knowing the 
H field, the respective components of E-field can be derived 
from the impedance condition: H = [z x E \/tjte, where 77 
is the transverse impedance of the mode. In MRI, H z com¬ 
ponent does not contribute to MR signal, but it can be im¬ 
portant to consider for mode excitations. In our effectively 
oversized waveguide with a high-permittivity dielectric in¬ 
sert, in principle, there are eight supported TW modes that 
can be excited and imaged as listed in Table I. The number 
of the modes in the fixed geometry setup can be further con¬ 
trolled by changing the permittivity of the dielectric inside 
the tube. In addition to hollow waveguide modes, the high 
reactance dielectric boundary conditions allow hybrid modes 
EH (or HE), which consist of linear combinations of TE and 
TM field components and therefore E Z ,H Z ^ 0. A com¬ 
bined system of a hollow metal guide that is partially filled 
with a concentric dielectric provides a radially anisotropic 


medium, and generally is not solved analytically 22 . Here, 
we consider a practical case of a quasi-uniform dielectric 
filling, where the dielectric-metal air gap is small compared 
to the wavelength and the Hi field is approximated by a 
field of a uniform waveguide. 

The rf modes could in principle be excited with a set of 
properly positioned elementary probes such as electric and 
magnetic dipoles. Due to wide availability of inductive loop- 
coils (magnetic dipoles) for near-field MRI we use these coils 
as probes to inductively couple into the dielectric guide. 
To predict the excitation modes these probes have to be 
inserted inside the dielectric, in practice however, we place 
them into the fields extrema of the specific mode outside of 
the dielectric . We consider two distinct cases: loop placed 
concentric (x-y plane) with the waveguide and orthogonal to 
the front edge of the waveguide (y-z plane). In the first case 
the magnetic field H||z couples into TE mode (Eq.(l)). In 
the second case Hlz and the probe predominantly excites 
lowest TEn mode. 

In general, the relationship between the MR image inten¬ 
sity and magnetic field map is non-trivial and specially de¬ 
signed MRI sequences exist to provide magnitude and phase 
magnetic field maps ■ . Such sequences usually target spa¬ 
tial profiles of H+ and separately, which are the right- 
and left-hand polarized components of Hi that correspond 
to transmit and receive sensitivity profiles of the rf coil re¬ 
spectively. For our application we are interested in obtain¬ 
ing actual propagating Hi field map from the MR image. 
It turned out that the MR signal intensity could provide a 
qualitative magnitude map of the Hi field for some impor¬ 
tant cases. For a gradient-recalled echo imaging sequence 
the signal intensity in each voxel is given by 

s cx psm(V\H+\^r)\(H^n (2) 

where we neglected relaxation effects; p is spin density dis¬ 
tribution, 7 is the gyromagnetic ratio, r is the rf excitation 
pulse duration, V is a dimensionless normalization factor 
proportional to coil driving voltage. Furthermore, in a small 
flip angle regime, which typically holds for TW MRI, the 
Eq. (2) can be rewritten as s oc \H+\\(H^n Next we use 
expressions for linearly driven coil 

H+ = l -{H lx + ifhy), H- = l -(H lx - iH ly y. (3) 

Substituting Eq. (3) into s we obtain 

8<x\\H* x + H* v \. (4) 


TABLE I. Propagating modes in the oversized circular waveg¬ 
uide with dielectric insert: mode index, critical wavelength A c , 
and propagating period A = 1/h. 


Mode 

TEn TMoi TE 2 i TE01, TMn TM21 TE12 TM02 

A c (cm) 
A(cm) 

15.4 11.8 9.3 7.4 5.5 5.3 5.1 

5 5.3 5.6 6.3 9.8 11.3 13.8 















3 


We further notice that the phase of the complex single mode 
TW field components (see Eq.(l)) is the spatial TW phase 
along z-axis, —ihz, so we can write Eq. (4) as 

s (X \\{H 2 1x + H 2 ly )e~ i2 ^~^\ = \\Ihl 2 , (5) 

where 0 is a fixed phase (ir/2 for TE modes), H\ x ,y are mag¬ 
nitudes of the respective field components with no z depen¬ 
dence. Notice, however, that the reflected wave gives rise to 
a standing wave pattern along z-axis and s oc cos(2hz). So, 
from the signal intensity distributions of the two orthogonal 
(/.e., coronal and axial) slices one can obtain a qualitative 
map of the magnitude of the magnetic field distribution. 
Eq. (5) is extended to a standing wave regime that includes 
a mixture of modes. However, Eq. (5) is not in general valid 
for the pure TW regime with multiple modes propagation 
and Hi MR mapping sequences are required to obtain the 
magnetic field distribution. 

III. EXPERIMENTAL RESULTS AND SIMULATIONS 

We performed proof of principle experiment on a 16.4 T 
horizontal-bore magnet interfaced to an Agilent DirectDrive 
console (Agilent Technologies, CA) with the resonance fre¬ 
quency of 698 MHz (A=43 cm). The magnet bore diameter 
is 26 cm, which is further reduced with a 12 cm gradient coil 
insert capable of reaching 1000 mT/min 150 fi s (Resonance 
Research, Inc., MA). The typical setup of a model waveg¬ 
uide with inductive coupling (Fig. 1) consists of a dielectric 
waveguide made of an acrylic tube (L=34 cm, outer diam¬ 
eter 10 cm, inner diameter 9 cm) and filled with deionized 
water or saline (£ r =80, a = 5 x 10 -4 S/m or a = 0.6 S/m 
respectively) and a simple transmit-receive loop-coil (D=7 
cm) placed at the edge of the tube. The respective modified 
wavelength in the dielectric is A e =4.9 cm. The orientation 
of the loop with respect to the edge of the tube (axial or 
longitudinal planes) defines the specific single or multiple 
rf modes excitations in the dielectric tube. The MR im¬ 
ages were obtained with a gradient-recalled echo sequence: 
FOV=20x20 cm, TR/TE=100/2 ms, flip angle set at 60°, 
10 (coronal and axial) slices with slice thickness of 5 mm, 
encoding matrix=128xl28, and bandwidth=89kHz. 

In our experimental study we used two excitation meth¬ 
ods by a single loop-coil that was oriented in two orthogonal 
planes: axial and coronal (longitudinal). Figure 2 shows MR 
images obtained for axial coil orientation in a sample of non- 
conductive dielectric (deionized (Di) water). Since the far 
edge of the dielectric tube represents unmatched boundary 
we expect a reflected TW thus creating an interference pat¬ 
tern in the FOV as shown in Fig. 2. For the actual magnetic 
field imaging purposes obtaining the interference fringes is 
complimentary for characterizations of the mode by mea¬ 
suring its period A. On the other hand, to eliminate back 
reflection and interference fringes we filled the tube with 
conductive dielectric (water) in otherwise same unmatched 
boundary setup. The resulted MR image is showing a sig¬ 
nal attenuated along z-axis (Fig. 3). In the latter setup the 
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FIG. 2. MR images (intensity in arbitrary units) of Di-water 
tube: axial (a), coronal off-axis (b) and on-axis (c). Dark 
spot in the center (a) and horizontal strip (c) is the null of 
TEoi mode. Vertical stripes (b, c) are interference fringes 
with period A/2 « 3 cm. (Curved bright spots on the edges 
(b, c) and ringing are artifacts from gradient edges.) 



FIG. 3. No interference effect: MR images (intensity in ar¬ 
bitrary units) of saline phantom, coronal off-axis (a) and on- 
axis (b). Saline conductivity strongly attenuates the reflected 
wave resulting in no interference fringes in this medium (the 
probe is located on the right side). Here, the gradient-related 
edge artifacts are amplified with respect to lower intensity 
decaying signal. 



FIG. 4. Simulations of Hi field of dominated mode TEoi in 
a partially filled dielectric waveguide: a) axial, b) coronal off- 
and c) on-axis slices, matching the intensity pattern in Fig. 2. 
Arrows indicate Hi vectors, the colorbar shows normalized 
magnitude of Hi. 

exact period of the mode could be derived from the signal 
phase plotted vs z coordinate. 

To validate obtained MR images and characterize mag¬ 
netic field patterns, we carried out simulations of the modes 
inside the screened dielectric guide using finite elements 
method (FEM) solver COMSOL Multiphysics (Burlington, 
MA). The results of the full-wave eigenmode Hi simulations 
of the same setup as used in Fig. 2 are shown in Fig. 4. 
The axial cross section view of the simulations shows ar¬ 
rows that depict direction of Hi and correspond to TEqi 
mode. Notice the extended view of the field pattern that 
covers otherwise MR invisible air all the way to the metal 
(gradient) boundaries. The simulations of Hi field profile 












4 



FIG. 5. Experimental MR images of Di-water tube with exci¬ 
tation of mixture of modes by the loop-coil oriented in coronal 
plane: coronal (a) and four axial slices (b). (Here, the inten¬ 
sity scale is arbitrary.) Simulations of Hi field of mixture of 
modes in a partially filled dielectric waveguide: c) coronal and 
d) four axial slices of mixture of modes (TMn and TEn)\ the 
colorbars show normalized magnitude of Hi. 


show excellent agreement with experimental MRI images. 

In another imaging example we used an alternative con¬ 
figuration of the loop coil that was positioned in orthogonal 
plane with Hi _L z and a dielectric tube filled with Di-water. 
The experimental results are shown in Fig. 5(a,b). Simi¬ 
larly to the first imaging example we obtained a longitudinal 
standing wave pattern (Fig. 5(a)) but with simultaneous ex¬ 
citation of a mixture of modes. The axial cross sections at 
various position along z (Fig. 5(b)) show different field pat¬ 
terns implying a mixture of at least two distinctive modes. 
The corresponding simulations of the Hi standing field pat¬ 
tern (Fig. 5(c,d)) matches the MR images and implies two 
basic modes: TEn and TMn. The nulls on z-axis (coronal 
slices) correspond to the standing waves of modes due to 
reflection from the far end of the dielectric guide. Over¬ 
all, the imaging experiment produced similar MR pattern to 
simulations of a relatively complicated field pattern. In the 
experimental images (coronal slices in Figs. 2,3, 5(a)) the 
high intensity star-like edge brightening and image distor¬ 
tions correspond to the artifacts due to spatially unencoded 
volume wrap-up of the dielectric guide that extends beyond 
the gradient allowed FOV and non-linear regions of the gra¬ 
dient coil respectively. In order to eliminate these artifacts 
one could utilize two compartment dielectric guide with MR 
non-visible dielectric (ceramic or D 2 0 ) and limit FOV within 
linear region of the gradient insert. 


IV. CONCLUSIONS 

In conclusion, we demonstrated the experimental method 
for direct imaging of magnetic rf field distribution in metal- 
dielectric waveguides based on TW MRI at ultra-high field. 
The magnetic field that is generated in MR active medium 
by the rf pulse facilitates imaging of the modes of cir¬ 
cular waveguide and opens up a practical way of non- 
perturbed imaging of the rf field distribution inside vari¬ 
ous hybrid waveguides, including waveguide mode convert¬ 
ers and transformers. We showed theoretically that TW 
excitation allows qualitative mapping of MR images to the 
magnetic rf field pattern and experimentally validated the 
technique by obtaining TW MR images of the magnetic field 
distribution of the rf modes of circular waveguide filled with 
deionized water in a 16.4 T small-bore MRI scanner and 
compared the MR images with numerical simulations. We 
found that modern systems with ultra-high magnetic fields 
(> 7 T) are well-suited for the rf field imaging. Our method 
could be the only direct way to visualize the rf field in the 
UHF band inside the waveguide structures apart from the 
numerical simulations. 
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